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Abstract 

Nonrelativistic Hamiltonians with large, even infinite, ground-state degener- 
acy are studied by connecting the degeneracy to the property of a Dirac operator. 
We then identify a special class of Hamiltonians, for which the full space of de- 
generate ground states in any spatial dimension can be exhibited explicitly The 
two-dimensional version of the latter coincides with the Pauli Hamiltonian, and 
recently-discussed models leading to higher-dimensional Landau levels are ob- 
tained as special cases of the higher-dimensional version of this Hamiltonian. 
But, in our framework, it is only the asymptotic behavior of the background 
'potential' that matters for the ground-state degeneracy. We work out in detail 
the ground states of the three-dimensional model in the presence of a uniform 
magnetic field and such potential. In the latter case one can see degenerate 
stacking of all 2d Landau levels along the magnetic field axis. 
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1 Introduction 



Quantum mechanical Hamiltonians admitting highly degenerate energy eigenstates 
are rather special, and they can thus serve as useful theoretical models to explain 
some novel properties exhibited by physical systems under certain circumstances. 
The most notable example is provided by the integer quantum Hall effect in two- 
dimensional (2d) electron gas, for which the quantized, and infinitely degenerate, 
energy eigenspace structure of the 2d Landau Hamiltonian [I] is largely responsible. 
The study of the 2d quantum Hall effect, including both integer and fractional ones, 
has been greatly benefitted by the elegant analytic properties of the Landau level wave 
functions (not only in Euclidean space, but also in curved or topologically different 
backgrounds [2]). In this 2d Landau Hamiltonian, the vector potential describing the 
background magnetic field is Abelian. 

Recently, it has been noticed by various authors that there also exist higher- 
dimensional, especially three-dimensional (3d), Hamiltonians the eigenstate structure 
of which exhibit Landau-level- like degeneracy [3H6]. A distinctive feature from the 
3d Euclidean-space Landau systems of Refs. [JHS] is that the related Hamiltonian 
typically contains spin-orbit coupling term E x p • a (for a radial electric field E for 
instance), which can be accommodated within the standard 'magnetic Hamiltonian' 
by having the spin-^ particle couple to a non-Abelian SU(2) gauge potential of the 
form A = Exit, One may expect that some nontrivial physical applications utilizing 
the latter Hamiltonians come out in the future. 

It would certainly be useful to know something on the Hamiltonian 'form' that 
can give rise to such highly degenerate eigenstate structure, especially for the ground 
state of the system. In the case of a nonrelativistic particle described by a spinor 
wave function in n-dimensional Euclidean (configuration) space, we may take the 
Hamiltonian to have the general form 

ft=2^(p-p + P(x,S)-p + Q(x,S)) , (1) 

where x = {x , • • • x n ) denote position coordinates, p = —iKV the momentum differ- 
ential operators, and S = {S^; 1 < i < j < n} represent related spin generators. (If 
the Hamiltonian involves functions Vi(x, S) which are linear in Sij, one may describe 
the related interaction by defining appropriate background gauge fields associated 
with the group SO(n); but in our work, we will not follow this line.) We are here 
particularly keen on the ground state degeneracy, and to facilitate our discussion it 
will be assumed that the function Q(x, S) in ([I]) has been chosen so that the lowest 
eigenvalue of % be equal to zero. Now our problem is : what kind of structures (and 
behaviors) for the functions "Pj(x, S) and Q(x, S) should be assumed to make our 
Hamiltonian to have a large, even infinite, zero-energy eigenspace? 

In the 2d case, Aharonov and Casher [7] made an interesting observation. If one 
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chooses 



Vi = -— Ai(x), (* = 1,2; x = (x,y)) 
c 

2 = (-V-A +? A.A)--. 3B (2) 

in JTJ, with arbitrary vector potential A(r) and its magnetic field B = d\Ai — d^A\ % 
one is led to the Pauli Hamiltonian 

W = (3) 

with 

H D = a • (p - -A) . (4) 

The it's are usual 2x2 Pauli matrices. Then, using the property of the 'Dirac 
Hamiltonian' Hp, they showed that ground state wave functions of H - the states 
with zero energy - can always be found analytically, with the degree of degeneracy 
determined by the total magnetic flux in accordance with the index theorem [8l 
[9]. Hence, if there exists a magnetic field approaching a nonzero constant value 
asymptotically, infinitely degenerate ground states follow always. 

In this work we will push this idea to higher dimensional setting to get a unified 
understanding on the problem posited above (in the Euclidean space only) . It is found 
that there exists a special family of Hamiltonians that enjoys the explicit ground-state 
integrability and encompasses all previously known models with the same property 
as special cases. New 3d models with interesting ground-state structure emerge along 
the way. In our 3d model with a uniform magnetic field, it is shown that entire 2d 
Landau-level wave functions, separated along the field direction, make up degenerate 
ground states. 

The rest of this work is organized as follows. In Sec. 2, a rather general discussion 
is offered on the structure of nonrelativistic Hamiltonians which can exhibit large 
ground-state degeneracy. In Sec. 3 we consider two classes of 3d Hamiltonians for 
which the full space of degenerate ground states can be exhibited explicitly as in the 
case of the 2d Pauli Hamiltonian. One is without any magnetic field and the other 
is with a magnetic field. In Sec. 4 we conclude with some remarks. In Appendix, we 
summarize the ideas in arbitrary spatial dimension. 



2 Hamiltonians with large ground-state degeneracy 

Our general strategy is as follows: large ground-state degeneracy results if a given 
nonrelativistic Hamiltonian % can be written in a square form ([3]), with the related 
Dirac Hamiltonian Hp having a structure that admits a large number of zero-energy 
eigenstates as required by the index analysis in open Euclidean spaces [9]. Here we 
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restrict to hermitian Hp so that the spectrum of the nonrelativistic Hamiltonian H 
is nonnegative. 

Being interested ultimately in nonrelativistic theory, we will write our Dirac 
Hamiltonian in the form 

H D = a • p + K(x, a,/3) , (5) 

by using hermitian Dirac matrices a = (7°7 1 , 7°7 2 , • ■ ■ 7°7 n ) an d /3 = 7° which 
satisfy the relations {a -1 , a 3 } = 2<5 y I and {/3,a*} = (with the matrices 7^ (fi = 
0, 1, • • • ,n) satisfying the Dirac-Clifford algebra relation {7^,7^} = —2rf LV ^ rf- v = 
diag(— 1, 1, 1, • • • )). The matrix (3 is taken to have the diagonal form 

'-(S-i)- < 6 » 

and the function K(x, ot,j3) in ©, a Dirac matrix polynomial, may be constrained 
by the condition {/3, K(x,a, f3)} = so that we have Hp} = and \P,"H] = 0. 
As P is diagonal, the Hamiltonian ([3]) can be written as 

With our choice of the Dirac Hamiltonian ([5]), eigenstates of the nonrelativistic 
Hamiltonian % of (J3j) can always be chosen to have a definite /3-parity, rj = ±1. For a 
given energy eigenstate Vl/ with the properties = and = where E > 0, 
another state ^ = Hpty will satisfy 

/ 3* = -ry§, K* = ^(^d) 3 * = M . (8) 

Hence all nonzero energy eigenstates of T~L appear pairwise, one with ry = +1 and 
the other with rj = — 1 ; this means that, when we leave out zero energy states, the 
two nonrelativistic Hamiltonians given in ([7]) are isospectral. For zero energy states, 
which must be the ground states, if exist, the situation can be different. 

Subsequent developments depend on the spatial dimension, and we shall study 
2d problem in this section and 3d problems in the next section. The case of higher- 
dimensional spaces will be considered in the context of specific models in the Ap- 
pendix. 

In 2d case, Dirac matrices are assumed by a 1 = Oi (i = 1, 2) and /3 = 03. Our 
Dirac Hamiltonian will then have the general form 

2 

H d ] = £ {*iPi + <n u i( x ) + mPV t (x)} . (9) 
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But we here have a 2d-specific identity ai/3 = —ie^aj, and so the t/j and Vi terms in 
([9]) are not independent. The form ([4]), used to define the 2d Pauli Hamiltonian, is 
essentially the unique possibility ; but, one may express this by an alternative form 

h d = Jy>i{K-*fciM*)) • (io) 

In the gauge V ■ A = 0, we can find a scalar 'potential' 4>(x) such that 

A x = --d24>, A 2 = -d 1( f> (11) 

e e 

and so the 2d magnetic field B = ^V 2 4>. The alternative form (|10|) becomes 

2 

#D } = ~ ih S ^ ^ + WW] • ( 12 ) 
i=l 

If we here take a background of the form (ft = -£-Bq{x 2 + y 2 ) as appropriate for 

(2) 

a uniform magnetic field -B(x) = i?o, the above Hamiltonian -f^ reduces to the 
2D version of a so-called Dirac oscillator Hamiltonian [TO] (but without mass term). 
(This was noticed also in Ref. [5]). If corresponds to a zero energy state of 
%V) = (H^) 2 /2M, it should satisfy the equation 



Hgho = -ihJ2°i ■ & + W) ^o = 0. (13) 

8=1 

With the P parity chosen such that f3^o = tj^/q (rj = ±1), we may set 

= e -^WF,(x) (14) 
to recast (I13p into the following equation on the 2-component spinor F^: 



Y t a i dl)F r , = 0. (15) 



i=i 

If we here write 

*Wi = ( f t ) ' F ^=- 1 = ( /°i ) ' (16) 

we then see from (fT5|) that the scalar functions /^(x) should satisfy 

(d x + ivd y )f v (x,y) = , (17) 

that is, /^(x, y) should be a function of the variable Q = x + irjy only. This was noted 
already in Ref. [7] , and one may use this result together with (fT4l) to produce ground 
state wave functions completely. 
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The number of independent ground states is decided by the asymptotic behavior 
of the scalar <j> which is in turn dictated by total magnetic flux on the 2d plane. If the 
magnetic field is asymptotically uniform, i.e., -B(x) — > Bq(> 0) as r = \J x 2 + y 2 — > 
co, 0(x) — > jf-B^r 2 asymptotically in the rotationally symmetric gauge. In this case, 
one finds infinitely degenerate ground states (all with the eigenvalue rj = +1), and 
the explicit ground state wave functions take the form 

{*,,(»)(*) = e-^(x + iy) n ; n = 0, 1,2, • • • } . (18) 

Note that all wave functions in (|18p remain normalizable as long as (lnr) _1 |</>(x)| — >• 
co ; with |</>(x)| — > clnr (c > 0), only a finite number of zero-energy modes are 
allowed. In the asymptotically uniform magnetic field, one could have chosen the 
asymptotically Landau gauge potential with </)(x) — > ^.B^y 2 as r — > co. Then, the 
ground state wave function would take the form of continuum states 

|* 0(fe) (x) = e -H*) e ik{x+iv) ; k real} . (19) 



3 Models based on 3d Dirac Hamiltonians 



In 3d case, Dirac matrices a 1 (i = 1,2,3) and f3 are provided by following 4x4 
matrices 

(20) 

Now our hermitian Dirac Hamiltonian Ho satisfying the condition {/3, Ho} = may 
have the form 



p + U + i/3V 



+ 7 5 W(x) + i 7 5 /3X(x) , (21) 



where 



7 5 = -zaW = ( J I 1 , (22) 



and all terms exhibited above are independent. Notice that our Dirac operator in 
(I2ip is entirely off-diagonal. Here the situation differs from the 2D case as Dirac 
spinors are defined in a 4-column space in contrast to nonrelativistic spinors defined 
in a 2-colunm space. Because of this, we are led to consider a pair of nonrelativistic 
Hamiltonians 1-l± (defined in 2-column spaces), according to 



Note that the subscripts ± denote the parity under 3d /3. In the form (I23j) we have 
natural candidate 3d Hamiltonians which can result in large ground state degeneracy 
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with suitably chosen background potentials L^(x), Vi(x), W(x) and V(x). Various 
potentials here may be considered in connection with their specific effects: J7j(x) for 
standard magnetic vector potentials, and Vi(x) for spin-orbit-coupling- like terms (see 
below), etc. 

In quantum-field-theoretic investigations the Dirac Hamiltonians similar to our 
form (f2Tj) have been discussed previously [TTHT3] ; there, potentials are typically those 
related to localized solitonic backgrounds (e.g., magnetic monopoles) and as such they 
usually involve internal space generators also. They also appear in some condensed 
matter literature discussing topological defects in insulators and superconductors; see 
Ref. [H] for instance. Zero modes of those Hamiltonians have been studied, together 
with responsible topological invariants. But, with the backgrounds of less restricted 
asymptotic behaviors (e.g. with unbounded potentials), and especially with an eye 
on their physical significance in nonrelativistic Hamiltonian contexts, the analysis is 
not so simple and up to our knowledge no systematic study has been made. Hence, 
leaving such to our future study, we shall below concentrate on some special class 
of Hamiltonians, which enjoy explicit ground-state integr ability and so can be used 
to exhibit some of the expected features. As we shall see, our model Hamiltonians, 
which have not been seriously considered in a relativistic setting, turn out to have 
some interesting nonrelativistic contents. 

Our first model in 3d is obtained from (|2ip by keeping only the potential V(x) 
in the form V(x) = — foV0(x), </>(x) being arbitrary. That is, we consider the Dirac 
Hamiltonian 

flg) = a ■ [p - ihpv</>(x)] , (24) 

and this model may in fact be considered in arbitrary spatial dimensions. In 3d case, 
a simple calculation using the Dirac matrices in (|20p yields the 3d Hamiltonian (|23|) 
which consists of a pair of nonrelativistic Hamiltonians 

= — (p 2 T 2h{V(f) x p) • <t + h 2 Vc/> • Vc/> T fr 2 V 2 c/>) • (25) 

If we here define the 'electric' field by E(x) = — V0(x), the second term in the 
right hand side of (|25l) obviously describes a spin-orbit coupling. Some recently 
discussed models with three-dimensional Landau-level-like structures correspond to 
special cases of this model, i.e., </>(x) = (const. ){x 2 + y 2 + z 2 ) in Refs. pH[5j, and 
0(x) = (const. )z 2 in the model of Ref. pj. In our discussion, however, the detailed 
profile of the potential 0(x) will be left largely arbitrary. As far as ground-state 
degeneracy structure is concerned, what matters is the asymptotic behavior of the 
potential - in our framework, a particular symmetry in the background potential is 
irrelevant. But our approach cannot say anything as regards possible excited-state 
degeneracy. 

The zero-energy eigenfunctions of the Dirac Hamiltonian (|24|) . as needed for 
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the ground states of the Hamiltonian (|25p . can be written as 

* (f) = e-^ x )F r/ (x) (26) 

with the P parity r\. Then the wave function F„(x) satisfies 

a ■ VF„(x) = . (27) 

For a given large- |x| behavior of </>(x), we usually obtain a physically acceptable wave 
function from (|27p only with a particular sign for r/. Aside from this, finding ground 
state wave function is now down to solving the background- free equation (127p . In 



this case, we again write F v=+ i = f ^ ~\ and F r)= _i = f ^ V but this time f v 

make 2-component spinors. Then the Dirac equation (|27p . with the Dirac matrices 
in (f20l) used, reduces to the 2-component spinor equation 



a • Vf„ = . (28) 

Since 

(o .. V ) 2 = ^ + a 2 + s 2 = (29) 

the 2-component spinor f should have harmonic functions as its components. Now, 
in view of the identity 

V 2 (z + ix cos u + iy sin it) p = , (ureal) (30) 

one may write the general solution to (I28p in the form |15| 

h(z + ix cos u + iy sin u, n)%(u)dn . (31) 

Here /i is an arbitrary function, and x( u ) represents a 2-spinor satisfying 

(o"3 + io"i cos n + ioi sin u)%(«) = , (32) 
whose solution has the following form 

X(«) = ( ) • (33) 



In pip any function /i is allowed as long as, when used in ()26|) . it gives rise to a 
physically acceptable ground-state wave function x &o( x ) • Given the potential with 
the asymptotic behavior </>(x) — > cr 2 (c > 0) as r = \J x 1 + y 2 + z 2 — > oo (as in the 
model of Refs. [US]), we must choose rj = +1 and then take for h following polynomial 
forms 

h = (z + ix cos u + iy sin u) l e imu , (I = 0, 1, 2, • • • ; m = -I, -I + 1, • • • I) . (34) 



7 



Then, using (|33p and performing the integration over u, we obtain the results for 
f(x) written using spherical harmonics : 



^±S±i r l Yi m (9,<j)) 
f + i(x)=( V ™ I (35) 



We thus have infinitely degenerate ground states in any such background. Actually 
all forms in ()35[) leads to normalizable wave functions as long as the asymptotics 
of the background potential is such that (lnr)~ 1 |g0(x)| — >• oo as r — > oo, and, if 
|<7</>(x)| — > clnr (c > 0) as r — >• oo, only a finite number of zero-energy modes survive. 
The detailed form of the background potential at finite r is not important for our 
discussion. 

On the other hand, with the background potential having the asymptotic behavior 
0(x) — y cz 2 (c > 0) as r — > oo (as in the model considered in Ref. [6]), we obtain 
bounded ground-state wave functions with /+(x) given by the continuum 

f+l(x) = e Hz+ixcosu+iy S inu) f \ ^ u real ) . ( 36 ) 

Note that, integrating this form over u with weight e irnu , one may also take for f + i 
the following expression: 



f +1 (x) = ( . «M ) ) , ( , real , ra = o, ±1, ±2, • • • ) (37) 

where x + iy = pe llp and J m {kp) denotes Bessel functions. 

For further discussions on these wave functions including possible physical ap- 
plications, see Ref. [3H6] . In 3d, the two expressions we have chosen above for the 
asymptotic form of the scalar <f> are not gauge equivalent as the scalar field is not 
connected to the gauge field in 3d and so the physics for these two cases are different. 

Let us now consider a little more complicated 3d Dirac Hamiltonian which has 
both gauge field and the scalar gradient as in 

H D = a-([p- - A(x)] - i^9V^(x)) (38) 

Under this choice, the related nonrelativistic Hamiltonians (1231) read as 



\ / P. P fxG \ 

H± = — ( [p--A] 2 T2W0x [p--A]-<7 <r-B + 7i 2 V0-V(^ft 2 V 2 0) . (39) 

2M V c c c / 

where B = V x A. If corresponds to a zero-energy ground state of H, it should 
satisfy the equation 

H d ^q = a-([-ihV- - A(x)] - ihpV^j ^0 = 0. (40) 
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With the P parity r\ chosen so that /3^o = V^o (v = il)> we may set 

* (f) = e -^WF,(x) (41) 
to recast (I40p into the following equation on F^(x): 

a ■ [ - ihV - ^A(x)] F„(x) = . (42) 

For a given large- |x| behavior of </>(x), we usually obtain a physically acceptable wave 
function from (|4ip only with a particular sign for r\. Aside from this, finding ground 
state wave function is now down to solving the simpler equation (|42p . Writing F + i = 
f+i\ ( 



^ y and F_i = ^ ^ J , ([12]) reduces to the 2-component spinor equations 

a ■ (V - -A)f„ = (43) 

For some magnetic flux, one could find normalizable zero-energy ground states if the 
scalar function cf> and the /3 parity r] were chosen suitably. 

Since a direct analysis of (|43p with an arbitrary 3d vector potential A(x) is 
impossible, let us restrict our attention to the case of a strictly 2d vector potential 

A = (A 1 (x,y),A 2 (x,y),A 3 = 0) (44) 

and rewrite as 

2 



(H%> + ip 3 )f(x, y, z) = 0, = ia 3 £ (Ti(pi - -A) (45) 

i=i c 

As ,Pz] = 0, we choose the eigenfunctions of p 3 as the solution of (f45j) . 



~ (2) 

Also note that H D is another hermitian Dirac Hamiltonian related to a 2d Pauli 
Hamiltonian, and so all normalizable energy engenfunctions of this operator can 
be found in principle. For the solutions of (I45p with (^3)' = we need zero-energy 

eigenfunctions of H D : they were found in Sec. 2. But, for the solutions with (p 3 )' 7^ 0, 

~ (2) 

we here need also nonzero-energy eigenfunctions of H D . As nonzero eigenvalues of 
this Dirac operator should appear pairwise with opposite signs, we may express the 
related eigenvalue equations schematically as 

H^h, = k e h e , H { * } a 3 h e = -k e a 3 h £ . (46) 

where £ labels all independent nonzero-energy engenstates. Then we can represent 
the (p 3 )' = ±ik£(^ 0) solutions to (|45p by the form 

f = e-^ z h e (x,y), f = e +! ^ z a 3 h e (x,y), (47) 
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Note that wave functions with imaginary eigenvalue of p% = —ihd^ would be out of 
consideration usually. However, in our case, there is additional factor e~ v ^ entering 
the wave function (|4ip . and an appropriate choice of r\ and (ft (imagine the scalar (ft 
having the asymptotic behavior (ft ~ cz 2 ) would provide sufficient fall off for finite ft£ 
at large \z\, making the ground wave function \&o m Q4ip normalizable. 

Especially, if the vector potential (|44p is chosen to be that of the uniform magnetic 



~ (2) 

field B = Bz, we can represent the eigenstates of H D in terms of the well-known 

Landau-level wave functions. Indeed, we then have H D (in the symmetric gauge) 
expressed as 



r(2) _ _ \ 1heB f a 
D ~ V c V a 

where a, a denote the creation and annihilation operators 



2|(a c - + a c -x). a = - ] /^(a ( -d (X ), (49) 

with ( = x + iy,Ai = —^eijdjX, and x = fUCC- The eigenstates and eigenvalues of 
are, respectively, 



, / \ / aW$ n \ , . l2heB 

Hn(x, y) = { _ ta \^ n j, h = ^ — (50) 

where £ = 0, ±1, ±2, • • • , and <3? n = Q n e~ x with n = 0, 1, 2, • • • , Especially, a choice 
(ft = cz 2 would make all 2d Landau levels be available for ground states of the 3d 
problem; but, as they are to be multiplied by the factor e ± ~h z ~ cz , the excited 2d 
Landau levels appear separated along the z-axis. As such degenerate stacking of all 
2d landau levels is possible, ground state degeneracy of this 3d system is infinite times 
larger than that of the 2d Landau system. The 2d Landau level functions in (|50p 
can also be used to find the ground state wave functions when the scalar field (ft has 
somewhat different asymptotic behaviors. In the case that we have (ft ~ T l-^l (t > 0) 
asymptotically, only a partial set of the above 2d Landau levels, that is, those with 
the Landau level index I satisfying the condition \kgjh\ < 7, would be acceptable for 
the ground states of the 3d system. 



4 Concluding Remarks 

In this paper we have shown that a fruitful way to study quantum systems exhibiting 
large ground-state degeneracy is to look for a connection to a Dirac operator. We 
are then led to a particular class of Hamiltonians, which exhibit explicit ground-state 
integrability (for any background potential with the given asymptotic behavior) and 
at the same time serve as a unified framework for some of the recently proposed 



10 



Hamiltonians in condensed matter physics. Needless to say, if the models considered 
in this work turn out to have some direct experimental relevance, it will be most 
welcome. To gain further insight, it should be desirable to have our analysis extended 
to the 3d model based on the full 3d Dirac Hamiltonian (|21|) . Considering our model 
in a curved space will be another interesting future problem. 
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A Arbitrary Dimension 

We shall comment on the case with a higher-dimensional-space version of the Hamil- 
tonian H obtained from the n-dimensional generalization of the Dirac Hamiltonian 
(124ft . For this, following Ref. [5], it is convenient to introduce rank-fc T-matrices, 1^ , 
% = 1, 2, • • • , 2k + 1 which are 2 k by 2 k and satisfy {Fi, Fj} = 25ij ; these matrices 
can be constructed iteratively using the recursive formulas 



r 



(AO 







-.(*-!) 



,(fc-l) 







(i < k-l) , F ( S 



o -a \ r ( fc ) / i o 

il J ' 2k+1 ~ \ -I 

(A.l) 

starting from rank-1 T-matrices T^ 1 "* = o"j (i = 1,2,3). Using these T-matrices, 
the SO(d) generators in the fundamental spinor representation can be given. In 
d = 2k + 1-dimensional space, one can take 



J 'j 



-S 



j 1 



T (k) p(fe) 



1 < i < j < 2k + 1) 



(A.2) 



for the rotation generators. On the other hand, in d = 2/c-dimensional space, we have 
two inequivalent set of SO(d) generators, {»%} and {Sy}, which can be identified 
with 



Set 1 



Sij 



Set 2 : S' 



-Sji 



.,Q> 



(fc-i) r (k-ih 



\Ff , for j = 



, for 1 < i < j < 2k - 1 
2k and 1 < i < 2k - 1 

, for 1 < i < j < 2k - 1 
= 2k and 1 < % < 2k - 1 



, (A.3) 
• (A-4) 
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We can also specify our Dirac matrices a 1 and /3 in even- or odd-dimensional space 
as follows : in 2k-dimensional space, take the 2 k x 2 k matrices 

a* = T\ k) (i = !,-■■ ,2k) , /3 = rf fc Vi , (A.5) 

and, in (2k + l)-dimensional space, use the expressions (which are 2 k+1 by 2 k+1 ) 

r (fc) \ ,. _ „, ^ „ / I 



« 4 = r (*) ^ (i = l,-" ,2fc + l) , 0= ( o _! ] • ( A-Ci) 



Inserting these Dirac matrices in the definition for 



Hd = a 



p - ih(3V(j)(x) (A.7) 



we come up with the Hamiltonian T~L = 2M^D in arbitrary spatial dimensions. In 
term of two components H± in ©, we have explicitly 

d = 2k 

2MH+ = p 2 - AhSijidifypj + /i 2 V0 • V(p - h 2 V 2 (t) 

2MU- = p 2 + AhS-jidifypj + fr 2 V0 • V0 + ft 2 V 2 0, (A.8) 

d = 2k + 1 

2M"H + = p 2 - AhSijidifypj + ?i 2 V0 • V0 - ft 2 V 2 

2MH_ = p 2 + mSij{di4>)pj + n 2 V0 • V</> + ft 2 V 2 . (A.9) 

For = 1 with (|A.8j) . setting S\2 = S[ 2 = 5 produces the 2D Pauli Hamiltonian 
with Ai = —e^djcf) ; taking k = 1 in ()A.9P and setting 5^ = \e^ k ak leads to our 
earlier 3D expression. 

For the ground-state wave functions one may solve the zero-energy Dirac equation 
Hd^o = by setting ^0 = e~^^F. Then, with the definite choice of the /3 parity 

f+i ' 



rj = ±1, we again obtain the equation (|27p for F„. If we here write F„ = i 







and F T?= _i = f ^ J, f±\ in d = 2/c-dimensional {d = 2k + 1-dimensional) space 
will have 2 k ~ 1 columns (2 k columns) and must satisfy the following equations: 

/2k-l \ 



d = 2k : f r f d i + irfckj Sj = , (A.10) 

/2fc+l \ 

d = 2k + 1 : ^ rf } d; fr, = (for both 77 = ±1) . (A.ll) 
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Note that f„, in both even- and odd-dimensional spaces, should have harmonic func- 
tions as its components. Hence, generalizing (13ip . we may express the column func- 
tion f„ in higher dimensional space by an integral 

f„ = / h(x d + iu • x)x (u) dQ d -2 (A.12) 

Ju£S d ~ 2 

where u • x = v^x 1 + • • • + u d ~ 1 x d ~ 1 , u denoting a vector which can take values on 
the sphere S d ~ 2 : (u 1 ) 2 + ■ • ■ + (n D_1 ) 2 = 1, and Xr)( u ) represents a column vector 
satisfying the condition 

/2fc-l \ 

d=2k : \Tu l T ! l k - 1) Jx v H = -VXr l (u), (A.13) 

d = 2k + 1 : iutr i k) + r 2fc+i^ ^(u) = (for both rj = ±1) (A.14) 

in even-and odd-dimensional spaces, respectively. Note that ()A.13p and (|A. 14|) may 
be written as eigenvector equations involving SO(D) generator matrices, viz., 

/2fc-l \ 1 

d = 2k : I u l S i>2k 1 x +1 (u) = -- x+1 ( u ) , 

/2fc-l \ 

( E J X-i(u) = ~-X-i(u) , (A.15) 

d = 2fc + l : ^E^ 5 ^+i^X ± i(u) = -^x ±1 (u) . (A.16) 

So one can use group theory to fix x v ( u ) here |16j . Then, depending on the asymptotic 
behaviors of the background potential 0(x), one may consider polynomial or/and 
exponential types for the function h in (1A.12|1 . just as in 2d and 3d cases treated 
earlier. This way, spinor wave functions analogous to the form (|35p (but now involving 
hyper-spherical harmonics [16]) or to the continuum expression (I36p can be obtained 
; using these in (|26|) will lead to normalizable (or at least bounded) ground-state wave 
functions only if the background potential 0(x) has 'right' asymptotic behaviors. For 
some of the explicit expressions regarding these wave functions, readers may consult 
Refs. [IHS] where they are discussed in the context of specially chosen background 
configurations. 
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